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9. )  Evaluation  of  Inertia  terms. 

It  was  found  that  the  calculated  Reynolds  number  for 
separation  is  16  times  to  small  compared  with  the  experi¬ 
mental  evidence.  This  discrepancy  was  attributed  to  the 
linearization  of  the  inertia  terms  in  Navier-Stokes  equations 
Therefore  an  estimation  of  the  inertia  terms  is  of  interest. 

In  (24)  the  first  term  represents  the  linearized  inertia 
forces,  the  second  and  third  one  on  the  left  side  the  second 
order  inertia  forces.  Por  an  estimation  the  surface  integral 
of  (24)  over  the  cross  section  of  the  gap  will  be  considered 

jU  ay'  df  ♦  f(\ ^  ay,'  -  y,'  a  yf  )  df  »  v  Jr  aa  y'  dp 

Regarding  (33)  one  sees  that  the  integrals  of  the  first  order 
inertia  term  and  of  the  viscosity  term  are  zero.  The  surface 
Integral  of  the  second  order  inertia  forces  will  be  transfor¬ 
med  by  Greens  theorem  into  a  boundary  integral: 

dF  -  f  ( y,  f£-  -  )  “s 


Only  the  outer  boundary  with  radius  r1  contributes  as  at  the 
inner  boundary  V* •  0  ,  yj  .  0  .  Introducing  (33)  and  the  direc¬ 
tion  of  n  equal  to  that  of  r  one  obtains  when  regarding  the 
boundary  conditions  (44,47) 


el  U*  5*  v 


The  following  values  were  computed  for  the  two  numerical 
examples 

K  -  3  •  1o5  ,  g*  -  196#6o595  U* 

K  -  1o4  ,  gT  -  1972,6o36  U* 

Introducing  the  experimental  values  r^  -  23  ma,  rQ  •  21  mm, 

e  »  1  mm  one  obtains  for  K  -  3  •  1o^  as  integral  value 
12 

-  35  U  •  The  second  order  inertia  terms  show  to  be  ne¬ 
gative  so  that  the  linearized  inertia  term  seems  to  be  too 
large.  If  this  is  so  a  too  small  Reynolds  number  is  calcu¬ 
lated.  This  can  be  shown  by  introducing  the  dimensionless 
radius  and  the  corresponding  stream  function  y  .  Then 
from  (22) 

Vr  *Vt  ~  n  *ay 

is  obtained.  An  average  value  of  the  Reynolds  number 


is  given  by 

ft  m  d F _ 

It  is  Indeed  too  small  if  the  average  value  of  the  dimen¬ 
sionless  inertia  forces  is  too  large. 

An  other  procedure  for  the  control  of  inertia  forces  is 
to  Introduce  in  advance  smaller  inertia  forces  than  in  the 
proceeding  calculation.  This  is  done  by  dividing  the  gap  in 
two  annum .  In  the  inner  one  inertia  foroes  are  neglected 
as  they  are  anyhow  zero  at  the  inner  boundary.  In  the  outer 
annuHus  the  linearized  inertia  forces  are  considered.  The 
width  of  the  two  annulti  was  assumed  as  equal.  The  radius  of 
the  interface  will  be  denoted  with  r’  and  the  corresponding 
variable  y  with  y'. 


The  solution  for  the  inner  annullus  was  already 
presented  in  chapter  (5)»  It  will  now  be  written 

V;  •  •  f  c,r  ♦  c,r‘"  ♦  c,r*  ♦  c,  r  In  r )  cos  y 

*  ef  b 4r  *  t^r*4  *  btrl  *  b*r  In  r  ]  sen  <p  }  (54 

Index  i  refers  to  the  inner  annulus,  index  a  to  the  outer 
annulus.  With  r  is  denoted  the  dimensionless  radius  r/rQ. 
The  boundary  conditions  are 

it  r‘>  1 ,  Vr  «  O  .  V*  •  0  (55 

at  r  ■  r',  y  -  y',  Vw  •  -  V*«‘  *  Vr»  (55 

For  the  outer  annulus  one  has  the  solution  (33)  and  the 
condition  (46)  for  the  outer  boundary. 

Furthermore  there  must  exist  a  steady  connection  of 
the  velocities  on  both  sides  of  the  interface.  This  gives 
the  additional  boundary  conditions  at  r  »  r',  y  ■  y'  resp. 

H  .  r’,  y  -  y';  *  V<n  (  W-*  *  (57 

Each  boundary  condition  gives  two  equations  one  for  the 
sine  and  one  for  the  cosine  term  of  (54).  Therefore  the 
boundary  conditions  lead  to  16  equations  by  which  the  8 
constants  c^,  Cg,  Cy  c^,  b^,  bg,  by  b^,  of  ( 54 )  and  the 
8  constants  aQ,  a^,  a2,  a y  mQ,  m ^ ,  mg,  m^  of  (3;:)  should 
be  determined.  However  two  of  these  equations  show  to  be 
linearily  dependent.  These  equations  are  VW  -  V*ri  at 
r  «  r',  y  -  y'  resp.  They  are  linear  combinations  of  (57). 

As  the  purpose  of  this  calculation  is  merely  to  analyse 
the  influence  of  too  large  inertia  terms  on  the  Reynolds 
number  for  separation  a  further  boundary  condition  may  be 


introduced  which  also  reduces  the  inertia  forces.  This 
is  the  condition  that  at  the  interface  the  second  radial 
derivatives  of  the  tangential  velocity  components  on  each 
side  of  the  interface  coincide.  This  gives  two  more  equations, 
one  for  the  sine  and  one  for  the  cosine  terms. 


V****  "  ( 58 

With  this  additional  condition  which  has  no  physical 
meaning  the  linear  dependency  is  removed.  It  was  already 
pointed  out  that  a  lack  of  physical  meaning  is  here  of  no 
importance  as  in  this  calculation  only  the  effect  of  smaller 
inertia  terms  than  in  the  previous  calculation  should  be 
investigated. 

Introducing  (54)  and  (33)  into  the  conditions  (55,  56, 

57#  58,  47)  one  obtains  the  following  boundary  conditions 

4.)  r  .  4 
c«  -  tj,  ♦  3  Cj  ♦ 
c,  *  cL  4  c, 


2.) 

v  .  r’  ,  y  •  y’ 

c,  - 

clTul  4  3  Cjt'1 

♦  c„  llnr' 

♦  4  )  *  f‘( y') 

b,  - 

b,r’*1  4  3  bjr  * 

*  b#(lnr' 

4  4)  -  9'<y') 

c.  * 

C4r’'*  4  c,v’*  * 

CHlnr'  • 

f  IV) 
y'*  1 

b,  ♦ 

b,r"»  4  b,r"  4 

Iwr'  - 

9  <*’> 

y’«  w 

2ctr-J  #  6ctr‘  4  c,r"  •  f(y') 

2b,r*J  4  bb,r’  4  bv r'~4  -  9" <y’) 

4  6c,  -  .  fm  ( y  ) 

-6  b,r’**  4  6  b,  -  b>r"1  •  9"  (y’> 


b4  •  b,  ♦  3b,  4  b,  •  0 


b,  ♦  b. 


.  Q 


(5* 


3  )  y  -  £0 


see  (47) 


The  condition  for  separation 
-£<1/*  u  )  -  0 

at  r  «  1  leads  with  (51a,  54)  to  the  expression 

2  U*  -  l{T  2ca  *  6c,  ♦  Ocos  f  ♦  [2b4  ♦  6  b,  ♦  bjit'w  J  •  0  ( 6 1 

Numerical  calculations  were  carried  through  for 
K  -  3  .  lo  ,  S0*  o,2,  y’  =  o,1,  rf  «  1,1 

With  these  quantities  the  16  constants  were  evaluated  1 iom 
(59#  6o,  47).  The  following  results  were  obtained. 


aQ  -  o,^7c9  V 
a1  -  -  2,43oo7  V 
&2  m  4 1,86o^1  V 

-  -  2o 1,64799  V 
»Q  •  -  0,52414  V 
m1  -  o,624b5  V 
m2  »  -  1,21o6l  V 

-  1l8,391o2  V 


C1  -  -  5,21248  .V 

c 2  *  2o,2o62o  V 

-  -14,^372  V 

c4  -  7o,39985  V 

b1  =  -  5, 1op97  V 

b^  »  -  8, o8644  V 

-  -13# 14616  V 
b4  -  -  42,46519  V  . 


Now  from  (6 1 )  the  excentricity  e  with  which  separa¬ 
tion  occurs  can  bs  determined.  Introducing  the  two  angles 
*f  ■  0°,  9o°  one  obtains 


5.45455  - 
0,26161, 
5.45455  - 
0,26951. 


t  •  2o,8499l  -  0 
•  -  5.45021  HI 

t  •  20,25889  -  0 
e  -  5.61477  an 


In  the  preceeding  calculation  with  linearized  Inertia 
terms  covering  the  whole  gap  it  was  found  e  *  1  for  the 
same  values  of  /,  and  K  -  5  •  1o^  or  Re  -  1,3  *  1o^  reap. 
The  experiments  had  shown  that  the  Reynolds  number  for 
separation  was  calculated  16  times  too  small.  Now  with 
reduced  inertia  terms  that  means  with  the  assumption  that 
only  in  the  outer  half  of  the  gap  inertia  forces  are  ao- 
ting  e  -  5,61477  on  is  found,  which  result  indeed  corres¬ 
ponds  more  closely  to  the  experiments.  A  direct  compari¬ 
son  is  not  possible  as  experiments  were  performed  only 
with  a  maximum  excentricity  e  ■  3,5  sa  corresponding  to 

l  -  0,167.  For  separation  the  experiments  gave  Re  - 
A  _  h 

2,1  •  1o  or  K  ■  4,8  *  1o  reap.  This  quantity  relating 
to  •  -  3,5  mm  must  Indeed  be  larger  than  for  e  -  5,6  mm. 
This  is  in  agreement  with  the  calculation  which  gave 
e  -  5.6  mm  for  K  -  3  •  10-*  or  Re  -  1,3  *  lo^  reap.  One 
sees  that  the  order  of  magnitude  of  experimental  and  cal¬ 
culated  Reynolds  numbers  now  agrees. 


Table  II 

Estimation  of  Inertia  terms. 

First  calculation 

K  -  3.1o5}  .  Re  -  1,3.  e  -  1,o13*2  ram;  £«  o,o48668 

Second  calculation 

K  -3.1  Re  -  1,3-  1o3,  e  -  5,45o21  mm;  £«  o,26l6l 

Experiment 

u  k 

K  -  i).  1o  ;  Re  -  2, 1 .  1o  ;  e  -  1  mm;  t  -  o,o48 

This  shows  clearly  that  the  linearization  causes  the 
discrepancy  of  calculated  and  measured  Reynolds  numbers 
for  separation.  By  omitting  partly  the  Inertia  terms  in 
the  second  calculation  Indeed  a  satisfactory  agreement 
with  experiments  is  obtained. 


1  ° . )  Stability  proof  of  Couette  flow  with  regard  to  per 


turbatlon  waves . 


In  the  preceeding  theoretical  investigation  it  was 
found  that  separation  of  the  Couette  flow  may  occur 
when  the  flow  is  bounded  by  excentric  cylinders  the 
outer  rotating  the  inner  at  rest.  The  close  agreement 
of  the  calculated  Reynolds  number  for  separation  and 
the  experimentally  determined  Reynolds  number  for  the 
first  occurance  of  perturbations  confirms  that  with  ex- 
centric  cylinders  the  generation  of  turbulence  is  a  se¬ 
paration  effect.  It  may  be  stated  that  the  same  effect  may 
be  due  to  vibrations,  as  separation  is  affected  by  Iner¬ 
tia  forces  independently  of  the  means  by  which  they  are 
created  i.e.  by  stationary  or  nonstationary  motion.  Those 
circumstances  may  explain  the  generation  of  turbulence 
in  earlier  experiments  *n  fact  the  review  of 

the  earlier  experimental  work  does  not  exclude  the  con¬ 
jecture  that  excentriclties  and  vibrations  could  have 
been  present.  This  could  mean  that  in  the  absence  of 
excentriclties  and  vibrations  a  Couette  flow  with  ro¬ 
tating  outer  cylinder  and  the  inner  cylinder  at  rest 
should  be  stable.  But  this  conclusion  would  not  be  com¬ 
plete  if  there  could  not  be  given  a  direct  proof  of  the 
stability  of  such  a  "pure"  Couette  flow.  By  eliminating 
excentriclties  and  vibrations  as  far  as  possible  sta¬ 
bility  was  found  to  the  highest  speeds  of  revolutions 
experimentally  obtainable  [3].  The  theoretical  treat¬ 
ment  applying  the  method  of  small  perturbation  waves 
on  the  complete  Navier  Stokes  equations  leads  to  Bessel 
functions  of  complex  order  as  elgenf unc tions  [5j.  As 
the  zeros  of  these  functions  are  net  known  the  eigen¬ 
values  cannot  be  determined  in  this  way.  Therefore  the 
solution  was  derived  as  a  series  expansion.  This  expan¬ 
sion  contains  few  numerical  errors  which  fortunately 
showed  to  be  of  subordinate  Influence  on  the  numerical 
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results.  In  eq.2o  of  [5]  the  coefficients  -  16.67} 
y$. 33  have  to  be  replaced  by  5*56}  22.22  and  in  eq.  21 

o 

of  [5]  the  i  terms  are  positive.  The  corrected  equa¬ 
tions  will  be  given  with  the  following  Rotations.  The 
fraction  of  the  width  d  of  the  gap  and  the  radius  r 
of  the  outer  cylinder  la  denoted  by  6, 


the  wave  length  of  the  perturbation  by^  and  the  number 
of  waves  along  the  circumference  by  k 


the  ratio  of  the  propagation  velocity  c  of  the  pertur¬ 
bation  and  the  velocity  V  of  the  outer  cylinder  by  , 


f  -  IP,* 'IF;  *  $ 

where  jr„  ft  are  the  real  and  imaginary  part  of  £  , 

the  Reynolds  number  by  R  , 


Vr 

m  ■*' 

V 


Finally  the  abbreviations 


q 

r 

s 


k  R  62 


k  R  (1  -  ?r) 

k*[t «•-  irhr] 


are  introduced. 


With  this  the  real  and  the  lmagninary  part  of  the  se 
ries  expansion  are 

33,3  (Hi)  4  [do, 5  ♦  11,3  k2  4  5,5  q  J  62 
4  [77.7  ♦  22,2  k2  4  5,5  q]&3 

♦  [  -  I9,79l5  -  43,05  k2  -  2,083  q  -  o,693  (k4  4  k2q) 

-  0,3472  (q2  -  p2)  ]  6  4 
+  [-  2,43o5  -  4,l5  k2  4  2,7  q  ♦  3,47?  k4 

4  2,o83  k2q  4  o,1157So  q2  -  0,578703  p2  -  0,4525  ps 


and 

5,5  +  22,2  6  -  (25,4525  4  2,33*5  k2)  52 
+  (11,34259  4  6,71255  k2)63 

-  (1,553  ♦  3,oo525  k2)64 

♦  (r f“  11**1  *  5,5  6  4  (9,72  4  1,35  k2)62 

-  (7,635  4  4,86335  k2)65  4  (2,7  4  2,083  k2)64] 

-  q62  [  1,553  -  2,083  6  4  o,65*  62 

4  (-  1,35  4  1,15755  6  -  o,23T55  62)J. 


There  exists  stability  when  §f;  <  0,  q<0  resp.  For  diffe 
rent  width  of  the  gap  and  different  wave  numbers  k  the 
following  stability  regions  for  jgr  were  calculated.  It 
has  to  be  mentioned  that  >1  has  no  physical  sense. 


Table  in 


Regions  of  in  which  <  0 


t 

K2 

min 

%r  max 

0,0l 

1 

o,51722 

1,28571 

loo 

0,51586 

19 

looo 

o,5o339 

It 

loooo 

o,36239 

It 

looooo 

0 

H 

0,1 

1 

o,64988 

1,29412 

loo 

o,55012 

it 

looo 

0 

it 

loooo 

0 

n 

0,2 

1 

0,76593 

1,25391 

loo 

o,47515 

R 

looo 

0 

It 

loooo 

0 

W 

o,4 

1 

°,9ol71 

1,1J31o 

loo 

0,58857 

looo 

0 

w 

loooo 

0 

It 

o,6 

1 

0,97159 

1,09525 

loo 

0,85994 

It 

looo 

o,8l647 

It 

loooo 

0,81026 

It 

oo 

o,8o954 

ti 

o 

« 

00 

1 

0,99654 

1 , 0^050 

loo 

o,977o5 

looo 

0,97399 

It 

loooo 

0,97363 

It 

90 

o,9735 9 

It 

0,9 

1 

o,99957 

1 , 00866 

loo 

0,99362 

It 

looo 

o,99295 

It 

loooo 

o,99288 

It 

oo 

o,99287 

It 

One  sees 

that  the  larger  the  width  of  the  gap 

the  smaller 

the  stabi 

lity  region  of  (fr 

which  in  the  limiting  case 

6  -1  seems  to  converge  to 

$r  -  1.  It  will  be 

shown  later 

that  this 

is  in  agreement  with  the  exact  solution,  which 

velocity  of  the  perturbance 


c  -  4  (59 

will  be  introduced.  As  mentioned  before  the  number  of  per¬ 
turbation  waves  around  the  circumference  is 

k  -  «  r  (60 


Introducing  the  perturbation  (53)  Into  Navier-Stokes 
equations  one  obtains  [11]  by  linearization 


W-  i&Ut  -  &  ]*»-*?■ ‘J  ♦  «-«'♦  £,  ♦ ,■&*) 
- ■  if  [(£/-*»>■ -  Hr  -*>■ - ^ 


(61 


_ £L  -ivl 

rf'-Jj*  rfr-y)  r  J 


This  expression  contains  the  complete  frictional  terras. 
An  exact  solution  of  eq.(6l)  will  be  deduced  . 

Introducing  (56)  one  sees  that  the  f»-term  on  the  left 
hand  side  vanishes.  Then  introducing  the  operator  [5] 


L(y)  * 


(62 


one  obtains 


(u-^c)L  - 


(63 


Thus  the  fourth  order  equation  (61)  is  reduced  to  the  two 
second  order  equations  (62,63). 

Introducing  the  dimensionless  quantities 


In  (62,63)  one  obtains 


(65 


L<V  -  t"~ 

jfe  +  f~o]L  *  0  (66 

Substituting 

tg*(/-H)*i*R(fS-i)  (67 

as  new  Independent  variable  (66)  is  transformed  to  the 
Bessel  equationj 

f 

tL'+lL  +  U*-*')!-  s  0  (68 

with  the  solution 

L  •  cJ/*)  +  ctH'<i)  (69 

(65)  i«  transformed  by  (67)  to 

*  0  (70 

Adding  (68,7o)  one  obtains 

A/.*  +  +  +)  m  o 

with  the  solution 


L  +  (f  -  <?,£*  +  e*z‘* 


Thus  tbs  complete  solution  of  (61)  is 


£'*-  c,  JKd)  -  eM  H*<*> 


(71 


with  constants  c^,c2»e^»c^. 

Ths  boundsrlss  srs  y  •  0,  y  •  r  or  with  (64,67) 

l*  )iKR(f-0  t  I,*0  (72 


Ths  boundary  Conditions  are 


<f*0>  £**0 


-k  i 

As  z  ,  are  Infinite  at  z  »  0  one  has 


c4  "  c2  *  © 


(73 


The  boundary  conditions  at  yield  the  homogenous 
equations 

t  =  't  l? - *  0 

7?  •  [$  «*,'**  3j  ’  0 

Putting  the  determinant  to  zero  one  obtains 


*-/ 


♦  *£ 


0 


(74 


This  expression  for  the  eigenvalues  z ^  can  be  transformed 
with  the  differential  formula  for  Bessel  functions 


£.  7  -7 

£  /w 


(75 


to 


As  the  Bessel  function  J  of  first  kind  has  only  real  roots 
the  eigenvalues  z ^  are  real.  This  requires  according  to 
(67) 

$r-  1  -  0  ,  0  (77 


where  are  the  real  and  Imaginary  parts  of  £  accor¬ 

ding  to  (64).  This  means  that  a  perturbation  once  origina¬ 
ted  can  only  rotate  with  the  angular  velocity  of  the  ro¬ 
tating  fluid.  £f-  <  0  means  damping.  Thus  a  flow  representing 
a  rigid  body  rotation  Is  stable.  This  confirms  the  former 
result  obtained  by  series  expansion  (s.  Tablelll).  It 
showed  that  the  larger  the  width  of  the  gap  the  more  the 
stability  region  converges  to  £r  -  1. 

The  solution  -  0  of  (76)  has  to  be  excluded  as 
2  «  0  represents  the  inner  boundary  according  to  (72). 

The  exaot  stability  proof  also  can  be  given  If  a  li¬ 
quid  annullus  of  the  width  d  is  rotating  as  a  rigid  body. 
The  boundary  is  the  same  as  In  the  proceeding  case, 
s.  (72)  but  z2  is  now  not  zero.  Therefore  the  constants 
c^.Cg  in  (71)  are  now  not  zero. 

The  coordinate  of  the  inner  boundary  io 

i2  -J-)  (7£ 


With 

'-f**  (79 

one  has 


(60 


o<  is  real. 


Introducing  (71)  into  the  two  boundary  conditions  (73) 
for  each  boundary  one  obtains 

c}  z*  +  C¥  z~*  -  c,  fait,)  -  Cj  Hk<2,)  =  0 

C j  K£,*~  -  cr4  C,  fyl,)  -  =  0 

c3  if  +  c4  zf*  -  c,  -  CjHxds)  -  0 

C3  Kzf~'~  C4  K  zf*~-  C,  -C2Hk(Z2)  c  0 

The  condition  that  the  determinant  must  be  zero  yields  if 
the  Wronski  determinant 

7k<*>  Hk<1>  -  -  |j-  (81 


is  introduceu 


-fcvK'aj .  .  j^zj  h^zj)] 


=  ( 


8K 

nz,i2 


By  rearranging  one  obtains 


With  the  differential  formulas  (75)  and 


(82 


which  hold  for  all  Winds  of  Bessel  functions  this  expres¬ 
sion  can  be  transformed  to 


2.  [~ 2-i)(7k( ^  *  z,  7k<£‘))  “  /W*>(7k  *  i2  "Jk 

+i?z,*[-HKH<t,H7*<M,> + 1 7Ktzi)  -  hI,u,>(7k<V  ~  fyM)] 


Here  again  introducing  (75# 32)  one  obtains 


*Z,  +  ^K-l ! s  0 

Expressing  the  Hank el  function  by  J  and  the  Neumann  func¬ 
tion  N  one  has 


Introducing  (80),  multiplying  the  whole  expression  by  a( 
and  differentiating  with  respect  to  •<  one  has 

*  0 


By  rearranging  and  introducing  (75*82)  one  obtains 


^ k<*z)[«*7k+iW  -  *  SNk-iU,)]  *  0 

As  z2#  0  one  has  the  condition 


JKH(z,) "  ^p^pA^v/Z/) 

7*  UjNk-iU,)  -NrJlJjK./U,) 


Aj  is  real  this  condition  means  that  the  ratios  of  the 
real  and  the  imaginary  parts  of  the  nominator  and  denomi¬ 
nator  must  be  equal.  This  is  not  possible  as  the  nomina¬ 
tor  and  the  denominator  are  of  different  degree  in  the 
lowest  power  of  z.  Therefore  the  condition  only  can  be  satisfil 
with  a  real  Independent  variable  z.  This  again  means  sta¬ 
bility.  It  can  be  seen  easily  that  (83)  reduces  to  (76)  if 
ci  tends  to  zero.  Then  the  JK  terms  vanish  and  ^(zg) 
can  be  cancelled.  What  is  left  is  (75). 


The  preceedlng  report  No.  1  of  June  196 1  was  devoted 
to  the  question  how  turbulence  may  be  generated  in  a 
Couette  flow  between  excentric  cylinders  when  the  outer 
one  is  rotating  and  the  inner  cylinder  is  at  rest. 

This  question  arose  when  recent  experiments  performed 
by  the  author  had  shown  definitely  stability  up  to  considerable 
speeds  of  revolution  f 5 J  whilst  the  earlier  experiments 
[l#2,3#*0  clearly  had  shown  transition.  As  in  the  new 
experiments  excentricities  and  vibrations  were  avoided 
as  far  as  possible  in  the  former  experiments  disregar¬ 
ded  excentricities  and  vibrations  could  have  affected 
the  transition  to  turbulence.  In  fact  a  review  of  the 
earlier  experimental  work  does  not  exclude  this  conjec¬ 
ture.  It  therefore  seemed  worthwhile  to  investigate  the 
generation  of  turbulence  by  excentricities  and  vibra¬ 
tions.  Turbulence  would  occur  here  as  consequence  of  se¬ 
paration.  As  separation  is  coupled  with  inertia  forces 
independent  of  the  means  by  which  they  are  produoed  it 
is  sufficient  to  consider  excentricity .  Restricting  to 
small  excentricities  with  regard  to  the  width  of  the 
gap  the  outer  boundary  which  is  regarded  as  excentric  can 
be  replaced  by  a  fictitious  centric  boundary  representing  the 
tniian  of  the  actual  boundary.  The  rotating  outer  cylinder 
produces  radial  velocities  at  the  fictitious  boundary* 

Therefore  nonhomogenous  boundary  conditions  exist  ex¬ 
cluding  an  eigenvalue  problem  and  necessitating  the  cal¬ 
culation  of  the  velocity  profiles.  Doing  this  it  was 
found  that  separation  can  occur. 

However  the  calculated  Reynolds  numbers  for  separation 
are  too  small  compared  with  experimental  observations. 

The  discrepancy  showed  to  be  too  large  as  to  be  explained 
by  the  earlier  ocouring  instability  at  the  inflection 
point  of  the  velocity  profile.  The  insufficiency  of  a 
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linear  approximation  seemed  to  be  more  likely.  In  fact 
the  mean  velooity  tends  to  zero  at  the  inner  boundary 
so  that  linearization  here  is  not  anymore  Justified.  Ob¬ 
viously  linearization  implies  to  large  inertia  forces. 
Therefore  in  this  report  a  new  calculation  was  performed 
with  inertia  forces’  which  should  be  smaller  than  the  ac¬ 
tual  inertia  forces.  This  was  done  by  neglecting  the 
inertia  forces  in  the  inner  half  of  the  annullus  and 
linearizing  the  inertia  forces  in  the  outer  half.  Now 
a  satisfactory  agreement  with  experiments  was  obtained. 

It  may  be  mentioned  that  linearized  boundary  calcu¬ 
lations  along  corrugated  walls  had  shown  a  sensitive 
influence  of  the  corrugation  on  transition  [7]  in 
agreement  with  the  result  obtained  here  that  the  calcu¬ 
lated  Reynolds  number  for  separation  i3  too  small.  Now 
as  a  comparison  with  experiments  is  available  this  sen¬ 
sitiveness  shows  to  be  more  attributed  to  the  insuffi¬ 
cient  linearization  than  to  physical  effects. 

Having  demonstrated  the  generation  of  turbulence  by 
an  excentricity  and  herewith  as  mentioned  before  also 
by  vibrations  a  straight  forward  proof  of  the  stability 
of  the  Couette  flow  with  rotating  outer  cylinder  on  the 
basis  of  propagating  perturbation  waves  is  still  of  in¬ 
terest.  The  exact  solution  could  be  derived  [5]  .  How¬ 
ever  the  eigenfunctions  are  Bessel  functions  of  complex 
order  which  are  not  yet  enough  explored.  Thus  a  series 
expansion  had  to  be  introduced  which  was  recalculated 
in  this  report.  It  shows  at  least  stability  regions. 
However  in  a  special  case  an  exact  stability  proof  can 
be  given.  It  is  the  rotation  of  a  fluid  as  rigid  body. 
This  is  also  the  limiting  case  of  the  centric  Couette 
flow  regarded  here  with  vanishing  radius  of  the  inner 
cylinder.  The  eigenfunctions  are  Bessel  functions  of 
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first  kind  the  zeros  of  which  are  the  eigenvalues.  As 
they  are  real  the  perturbation  waves  are  damped  and  ro¬ 
tate  with  the  body.  Damping  means  stability.  The  exact 
stability  proof  is  also  given  for  a  liquid  annullus  ro¬ 
tating  as  a  rigid  body.  Also  in  this  case  the  eigenva¬ 
lues  turn  out  to  be  real  so  that  stability  exists. 
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